Integrating microbial abundance time series with fermentation dynamics of the rumen microbiome via mathematical modelling.

The rumen represents a dynamic microbial ecosystem where fermentation metabolites and microbial concentrations change over time in response to dietary changes. The integration of microbial genomic knowledge and dynamic modelling can enhance our system-level understanding of rumen ecosystem's function. However, such an integration between dynamic models and rumen microbiota data is lacking. The objective of this work was to integrate rumen microbiota time series determined by 16S rRNA gene amplicon sequencing into a dynamic modelling framework to link microbial data to the dynamics of the volatile fatty acids (VFA) production during fermentation. For that, we used the theory of state observers to develop a model that estimates the dynamics of VFA from the data of microbial functional proxies associated with the specific production of each VFA. We determined the microbial proxies using CowPi to infer the functional potential of the rumen microbiota and extrapolate their functional modules from KEGG (Kyoto Encyclopedia of Genes and Genomes). The approach was challenged using data from an in vitro RUSITEC experiment and from an in vivo experiment with four cows. The model performance was evaluated by the coefficient of variation of the root mean square error (CRMSE). For the in vitro case study, the mean CVRMSE were 9.8% for acetate, 14% for butyrate and 14.5% for propionate. For the in vivo case study, the mean CVRMSE were 16.4% for acetate, 15.8% for butyrate and 19.8% for propionate. The mean CVRMSE for the VFA molar fractions were 3.1% for acetate, 3.8% for butyrate and 8.9% for propionate. Ours results show the promising application of state observers integrated with microbiota time series data for predicting rumen microbial metabolism.

The effects of residual energy intake on nutrient use, methane emissions and microbial composition in dairy cows.

For sustainable food production selection and breeding of feed efficient animals is crucial. The objective of this study was to evaluate whether multiparous dairy cows, ranked during their first lactation based on residual energy intake (REI) as efficient (low; L-REI) or inefficient (high; H-REI), differ in terms of nutrient use efficiency, methane emissions, rumen fermentation, and gut microbiota composition. Six L-REI and 6 H-REI cows were offered two diets with either a low or high proportion of concentrates (30 vs. 50% of DM) on two consecutive periods of 21 d. Gas exchanges, milk yield, feces and urine excretions were measured in open-circuit respiratory chambers. The results indicated that L-REI cows had higher methane yields (22.6 vs. 20.4 g/kg DM intake) and derived more energy (energy balance - 36.6 vs. - 16.9 MJ/d) and protein (N balance - 6.6 vs. 18.8 g/d) from the tissues to support similar milk yields compared to H-REI cows. Nutrient intake and digestibility were not affected by REI, and there were no interactions between REI and diet. Milk yield, milk production efficiency, and milk composition were not affected by REI except for milk urea concentration that was higher for L-REI cows (14.1 vs. 10.8 mg/100 ml). The rumen and fecal microbiota community structure and function were associated with both the diet and REI, but the diet effect was more pronounced. The current study identified several physiological mechanisms underlying the differences between high and low REI cows, but further studies are needed to distinguish the quantitative role of each mechanism.

Effects of dietary forage-to-concentrate ratio and forage type on milk phospholipids and fatty acid composition of polar lipids.

The effects of grass silage and red clover silage on milk fatty acid (FA) composition are extensively studied, but little is known of their effects on minor lipid constituents of milk fat globule membrane. We investigated the effects of forage:concentrate (FC) ratio in grass silage-based diets and forage type (grass silage vs. red clover silage) on selected molecular species of milk phospholipids (PL) and the FA composition of PL. Ten multiparous Nordic Red cows were offered following dietary treatments: grass silage-based diets containing 70:30 (HG) or 30:70 (LG) FC ratio or a red clover silage-based diet (RC) comprising 50:50 FC ratio on a dry matter basis. The most abundant molecular species within the phosphatidylcholines was 16:0-18:1 phosphatidylcholine that was increased by 18% in HG compared with LG milk. Dietary treatments did not affect the relative proportion of 18:1-18:1+18:0-18:2 phosphatidylethanolamine that was the most prevalent species (ca. 44%-45%) in that class. We identified the d18:1-22:0 sphingomyelin as the most abundant sphingomyelin species that tended to increase in HG milk compared with LG. The FC ratio did not affect the relative proportions of saturated FA nor monounsaturated FA in PL, but the proportion of cis-9 18:1 was elevated in HG versus LG milk, whereas the proportion of 18:2n-6 was 50% higher in LG versus HG milk. The RC diet increased monounsaturated FA and 18:3n-3 levels in PL compared with grass silage-based diets and decreased the relative proportion of saturated FA. However, the RC diet did not affect the relative proportion of polyunsaturated FA in PL, although red clover silage typically increases the proportion of polyunsaturated FA in milk fat. This study provides valuable knowledge of the minor lipid components in milk on species level in relation to common feeding strategies in high-forage systems.

Long-term effects of early-life rumen microbiota modulation on dairy cow production performance and methane emissions.

Rumen microbiota modulation during the pre-weaning period has been suggested as means to affect animal performance later in life. In this follow-up study, we examined the post-weaning rumen microbiota development differences in monozygotic twin-heifers that were inoculated (T-group) or not inoculated (C-group) (n = 4 each) with fresh adult rumen liquid during their pre-weaning period. We also assessed the treatment effect on production parameters and methane emissions of cows during their 1st lactation period. The rumen microbiota was determined by the 16S rRNA gene, 18S rRNA gene, and ITS1 amplicon sequencing. Animal weight gain and rumen fermentation parameters were monitored from 2 to 12 months of age. The weight gain was not affected by treatment, but butyrate proportion was higher in T-group in month 3 (p = 0.04). Apart from archaea (p = 0.084), the richness of bacteria (p < 0.0001) and ciliate protozoa increased until month 7 (p = 0.004) and anaerobic fungi until month 11 (p = 0.005). The microbiota structure, measured as Bray-Curtis distances, continued to develop until months 3, 6, 7, and 10, in archaea, ciliate protozoa, bacteria, and anaerobic fungi, respectively (for all: p = 0.001). Treatment or age × treatment interaction had a significant (p < 0.05) effect on 18 bacterial, 2 archaeal, and 6 ciliate protozoan taxonomic groups, with differences occurring mostly before month 4 in bacteria, and month 3 in archaea and ciliate protozoa. Treatment stimulated earlier maturation of prokaryote community in T-group before month 4 and earlier maturation of ciliate protozoa at month 2 (Random Forest: 0.75 month for bacteria and 1.5 month for protozoa). No treatment effect on the maturity of anaerobic fungi was observed. The milk production and quality, feed efficiency, and methane emissions were monitored during cow's 1st lactation. The T-group had lower variation in energy-corrected milk yield (p < 0.001), tended to differ in pattern of residual energy intake over time (p = 0.069), and had numerically lower somatic cell count throughout their 1st lactation period (p = 0.081), but no differences between the groups in methane emissions (g/d, g/kg DMI, or g/kg milk) were observed. Our results demonstrated that the orally administered microbial inoculant induced transient changes in early rumen microbiome maturation. In addition, the treatment may influence the later production performance, although the mechanisms that mediate these effects need to be further explored.

Review: Problems in determining metabolisable protein value of dairy cow diets and the impact on protein feeding.

Accurate determination of dairy cow's metabolisable protein (MP) requirements is essential for maximising milk yield and minimising N input in dairy production systems. The main objectives of feed evaluation systems are to predict accurately the nutritive value of feed resources and production responses to ingested nutrients. In recent years, protein evaluation systems have made progress in complexity but our analysis indicated that in comparison to simple MP models, or even models based on CP and metabolisable energy intake, new models have failed to improve milk protein yield predictions. Previous meta-analyses have indicated low marginal efficiencies for incremental CP intake even with high-quality feeds as protein supplements. Most current protein systems tend to overestimate MP supply from rumen undegradable protein and overlook the effects on rumen microbial protein synthesis. These findings suggest that there is scope for improving predictions of MP supply. This review highlights and discusses some evident problems in the current MP systems that may partly explain the modest progress achieved over the last few decades.

Fresh Rumen Liquid Inoculant Enhances the Rumen Microbial Community Establishment in Pre-weaned Dairy Calves.

The development of the functional rumen in calves involves a complex interplay between the host and host-related microbiome. Attempts to modulate rumen microbial community establishment may therefore have an impact on weaning success, calf health, and animal performance later in life. In this experiment, we aimed to elucidate how rumen liquid inoculum from an adult cow, provided to calves during the pre-weaning period, influences the establishment of rumen bacterial, archaeal, fungal, and ciliate protozoan communities in monozygotic twin calves (n = 6 pairs). The calves were divided into treatment (T-group) and control (C-group) groups, where the T-group received fresh rumen liquid as an oral inoculum during a 2-8-week period. The C-group was not inoculated. The rumen microbial community composition was determined using bacterial and archaeal 16S ribosomal RNA (rRNA) gene, protozoal 18S rRNA gene, and fungal ITS1 region amplicon sequencing. Animal weight gain and feed intake were monitored throughout the experiment. The T-group tended to have a higher concentrate intake (Treatment: p < 0.08) and had a significantly higher weekly weight gain (Treatment: p < 0.05), but no significant difference in volatile fatty acid concentrations between the groups was observed. In the T-group, the inoculum stimulated the earlier establishment of mature rumen-related bacterial taxa, affecting significant differences between the groups until 6 weeks of age. The inoculum also increased the archaeal operational taxonomic unit (OTU) diversity (Treatment: p < 0.05) but did not affect the archaeal quantity. Archaeal communities differed significantly between groups until week 4 (p = 0.02). Due to the inoculum, ciliate protozoa were detected in the T-group in week 2, while the C-group remained defaunated until 6 weeks of age. In week 8, Eremoplastron dilobum was the dominant ciliate protozoa in the C-group and Isotricha sp. in the T-group, respectively. The Shannon diversity of rumen anaerobic fungi reduced with age (Week: p < 0.01), and community establishment was influenced by a change of diet and potential interaction with other rumen microorganisms. Our results indicate that an adult cow rumen liquid inoculum enhanced the maturation of bacterial and archaeal communities in pre-weaning calves' rumen, whereas its effect on eukaryotic communities was less clear and requires further investigation.

The effect of dietary forage to concentrate ratio and forage type on milk fatty acid composition and milk fat globule size of lactating cows.

We examined the effects of 2 grass silage-based diets differing in forage:concentrate (FC) ratio and those of a red clover silage-based diet on intake, milk production, ruminal fatty acid (FA) biohydrogenation, milk FA composition, and milk fat globule (MFG) size distribution. Ten multiparous Nordic Red cows received the following treatments: grass silage-based diets containing high (70:30, HG) or low (30:70, LG) FC ratio or a red clover silage-based diet with an FC ratio of 50:50 (RC) on a dry matter basis. Determinations of MFG were performed from fresh milk samples without addition of EDTA so the results of fat globules >1 µm in diameter are emphasized instead of the entire globule population. Lower FC ratio in grass silage-based diets increased milk production with no effect on daily fat yield, leading to 13% lower milk fat concentration. The effect of FC ratio on MFG size was moderate. It did not affect the volume-weighted diameter in grass silage-based diets, although LG lowered the volume-surface diameter of MFG in the size class >1 µm compared with HG. Compared with HG, feeding LG moderately decreased the biohydrogenation of 18:2n-6, leading to a higher level of polyunsaturated fatty acids in milk fat. Feeding RC lowered milk fat concentration and daily milk fat yield compared with grass silage-based diets. The volume-weighted diameter of MFG in the size class >1 µm was smaller in RC milk compared with grass silage-based diets. Feeding RC increased the flow of 18:3n-3 at the omasum by 2.4-fold and decreased the apparent ruminal 18:3n-3 biohydrogenation compared with grass silage-based diets despite similar intake of 18:3n-3. It also resulted in the lowest amount of saturated FA and the highest amounts of cis-9 18:1, 18:3n-3, and polyunsaturated FA in milk. In conclusion, LG decreased milk fat content and induced minor changes in MFG size distribution compared with HG, whereas RC lowered milk fat production, altered milk FA composition to nutritionally more beneficial direction, and led to smaller MFG compared with grass silage-based diets.

Dietary fish oil supplements modify ruminal biohydrogenation, alter the flow of fatty acids at the omasum, and induce changes in the ruminal Butyrivibrio population in lactating cows.

Four lactating cows fitted with ruminal cannulae and fed a grass silage-based diet were used in a 4 × 4 Latin square with 28-d periods to investigate the effects of incremental dietary fish oil (FO) supplementation (0, 75, 150, or 300 g/d) on the flow of fatty acids at the omasum and populations of rumen bacteria capable of biohydrogenation. FO decreased silage intake and ruminal volatile fatty acid concentrations and promoted an increase in molar butyrate and propionate proportions at the expense of acetate. Extensive ruminal biohydrogenation of 20:5(n-3) and 22:6(n-3) resulted in corresponding increases in numerous 20- and 22-carbon unsaturated fatty acids at the omasum. Omasal flow of several 20-, 21-, and 22-carbon all-cis (n-3) PUFA exceeded the intake from FO. Supplements of FO also induced a dose-dependent decrease in 18:0 and increased trans 18:1 and trans 18:2 flow at the omasum. Trans-11 was the major 18:1 intermediate in digesta, while FO induced quadratic increases in trans-10 18:1 flow, reaching a maximum of 300 g/d. FO had no substantial influence on omasal flow of CLA. Results suggest that one or more fatty acids in FO inhibit the reduction of trans-18:1 and trans-18:2 intermediates by ruminal microorganisms. qPCR based on 16S rRNA genes in omasal digesta indicated that key Butyrivibrio spp. declined linearly in response to FO. Dose-dependent increases in ruminal outflow of biohydrogenation intermediates containing one or more trans double bonds in response to FO has major implications for host metabolism and the nutritional quality of ruminant foods.

Ruminal infusions of cobalt-EDTA reduce mammary delta9-desaturase index and alter milk fatty acid composition in lactating cows.

Ruminal administration of a triple indigestible marker system comprised of cobalt EDTA (CoEDTA), ytterbium acetate (YbAc), and chromium-mordanted straw (CrS) decreases product:substrate ratios for Delta9-desaturase in bovine milk fat. This experiment was designed to identify the marker(s) responsible and develop an alternative system for simultaneous determination of nutrient flow in the gastro-intestinal tract and milk fatty acid composition. Five lactating dairy cows were used in a 5 x 5 Latin square with 21-d periods to evaluate the effects of YbAc, CoEDTA, and CrS independently or as part of a triple marker system (TMS), and CrEDTA as an alternative to CoEDTA on milk fat composition. Markers were administered in the rumen over a 7-d interval and samples of milk were collected on d -1, 3, 7, and 11. Both TMS and CoEDTA alone reduced the concentrations of milk fatty acids containing a cis-9 double bond, whereas YbAc, CrS, and CrEDTA had no effect. Reductions in product:substrate ratios for Delta9-desaturase were time dependent and evident within 3 d of administration. Ruminal infusion of CoEDTA for 7 d induced mean decreases in milk cis-9 14:1/14:0, cis-9 16:1/16:0, cis-9 18:1/18:0, and cis-9, trans-11 conjugated linoleic acid/trans-11 18:1 concentration ratios of 47.7, 26.7, 40.3, and 42.6%, respectively. In conclusion, ruminal infusion of CoEDTA alters milk fatty acid composition and appears to inhibit Delta9-desaturase activity in the bovine mammary gland. Results indicate that a TMS based on CrEDTA, YbAc, and indigestible neutral detergent fiber can be used for estimating nutrient flow without altering milk fat composition in lactating cows.

Transfer of absorbed cis-9, trans-11 conjugated linoleic acid into milk is biologically more efficient than endogenous synthesis from absorbed vaccenic acid in lactating cows.

Cis-9, trans-11, the major isomer of conjugated linoleic acid (CLA) in bovine milk fat, is derived from ruminal biohydrogenation of 18:2 (n-6) and endogenous conversion of trans-11 18:1 (vaccenic acid; VA) in the mammary gland. Most evidence to date suggests that endogenous synthesis is the major source of cis-9, trans-11 CLA, but the extent of VA desaturation is less well defined. Four lactating cows were used in consecutive 4 x 4 Latin squares to examine changes in milk fatty acid composition and secretion in response to abomasal infusions of lipid supplements enriched with cis-9, trans-11 CLA (88.8%) or VA (29.4%). Treatments were infused over 4-d, followed by a 3-d washout, during 7 d experimental periods and administered to deliver 0, 3, 6, and 12 g cis-9, trans-11 CLA/d (Expt. 1) or 0, 7.5, 15 and 30 g VA/d (Expt. 2). Infusions of cis-9, trans-11 CLA increased linearly milk cis-9, trans-11 CLA concentrations from 0.68 to 1.46 g/100 g fatty acids. Abomasal infusions of VA increased linearly milk VA and cis-9, trans-11 CLA content from 1.22 to 2.72 and 0.61 to 1.24 g/100 g fatty acids, respectively. Changes in milk fatty acid secretion indicated that 28.9% of VA was converted to cis-9, trans-11 CLA. Results provide evidence that conversion by Delta9-desaturase to cis-9, trans-11 CLA in the lactating cow is independent of postruminal VA supply. In conclusion, endogenous synthesis via VA was equivalent to approximately 21% of the response to increases in cis-9, trans-11 CLA available for absorption.

Determination of digesta flow entering the omasal canal of dairy cows using different marker systems.

Four studies were conducted to compare the effect of four indigestible markers (LiCoEDTA, Yb-acetate, Cr-mordanted straw and indigestible neutral-detergent fibre (INDF)) and three marker systems on the flow of digesta entering the omasal canal of lactating dairy cows. Samples of digesta aspirated from the omasal canal were pooled and separated using filtration and high-speed centrifugation into three fractions defined as the liquid phase, small particulate and large particulate matter. Co was primarily associated with the liquid phase, Yb was concentrated in small particulate matter, whilst Cr and INDF were associated with large particles. Digesta flow was calculated based on single markers or using the reconstitution system based on combinations of two (Co + Yb, Co + Cr and Co + INDF) or three markers (Co + Yb + Cr and Co + Yb + INDF). Use of single markers resulted in large differences between estimates of organic matter (OM) flow entering the omasal canal suggesting that samples were not representative of true digesta. Digesta appeared to consist of at least three phases that tended to separate during sampling. OM was concentrated in particulate matter, whilst the liquid phase consisted mainly of volatile fatty acids and inorganic matter. Yb was intimately associated with nitrogenous compounds, whereas Cr and INDF were concentrated in fibrous material. Current data indicated that marker systems based on Yb in combination with Cr or INDF are required for the accurate determination of OM, N and neutral-detergent fibre flow. In cases where the flow of water-soluble nutrients entering the omasal canal is also required, the marker system should also include Co.